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ABSTRACT 


The general relativistic equations of stellar structure and 
evolution are re formula ted in a notation which makes easy contact 
with Newtonian theory. Also, a general relativistic version of 
the mixing- length formalism for convection is presented. Finally, 
it is argued that in previous work on spherical systems general 
relativity theorists have identified the wrong quantity as "total 
mass-energy inside radius r." 


Subject headings: interiors: stellar, relativity. 


i 


I. INTRODUCTION 


The general relativistic equations of stellar structure for zero- 
temperature stars (neutron stars) were first presented in their modern form 
by Oppenhelmer and Volkoff (1939). Two decades later the discovery of 
discrete, galactic X-ray sources (Clacconi et al . 1962, Gursky et al. 1933) 
motivated theoretical studies of hot, relativistic neutron stars by Chiu 
and Salpeter (1964), Morton (1964), and Tsuruta (l964)j and the huge energy 
requirements of strong radio sources motivated Hoyle and Fowler ( 1963a, b) 
to develop the theory of hot, “fficiently convective supermasslve stars in 
which, it was soon realized, relativistic effects can be important (Feynman 
1964, Chandrasekhar 1964, Fowler 1964). In response to these developments, 
and others, Bardeen (1965), Misner and Sharp (1965), and Lindquist (1966) 
developed the theory of diffusive heat transfer in relativistic stars, and 
Bondi (19o4), Chandrasekhar (1965) and Thorne (1966a) elucidated the rela- 
tivistic version of the Schwarzschild criterion for convection. All of these 
pieces of relativistic stellar theory were put together and combined with 
relativistic equations for nuclear energy generation by Hameen-Anttila and 
Anttila (1966) and by Thorne ( 1966b, 1967) to give the currently standard 
version of the relativistic equations of stellar structure and evolution. 

Recently Anna &ytkow and 1 began analyzing the structure of red super- 
giant stars with degenerate neutron cores (see the following paper and 
references cited therein). For this purpose the standard relativistic 
stellar equations are unsatisfactory in two ways; (i) they do not make easy 
contact with the standard Newtonian equations; and (ii) they do not include 
a mixing-length formalism for convective energy transport. The purpose of 
this paper is to remedy these defects by (i) translating the relativistic 
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equations Into a new notation, and (ii) presenting a straightforward rela- 
tivistic generalization of the standard Newtonian mixing- length equations. 

No detailed derivations will be given because the translation from the old 

0 

notation (Thorne 1966b, 1967) to the new is straightforward; and the deriva- 
tion of the relativistic mixing-length theory is identical to the Newtonian 
derivation, if one works in the proper reference frame of a static relativistic 
observer. 

Throughout the paper c.g.s. units will be used; the speed of light c 
and Newton's gravitation constant G will not be set equal to unity. 


II. FUNDAMENTAL VARIABLES 


As our independent thermodynamic variables we choose the following — 
all of which are determined by measurements using standard, physical rods 
and clocks, in the mean local rest frame of the baryons. After the symbol 
■for each quantity, we Indicate its units in brackets. 


r_i^l ^ (density of "rest mass") ^ hydro jen \ 

Lcm^-I ^atom in its ground state J 

/number density\ 


y of baryons j * 


T[°K] 


s (temperature); 




fractional abundance of nuclear species i, \ 
by rest mass or equivalently by baryon numbery 


(la) 

(lb) 

(l c) 


As our independent radial and time variables we choose 




/"rest mass" inside\ _ / mass of one hydrogen \ 
^ "radius r" J ~ ^atom in its ground state] 


/total number of baryons 


y I — aryons'\ , 

y inside radius r j * 


(2a) 


o 


1 


, , /"Schwarzschild time\ / time coordinate such that Ib\ 

t[secj 5 ? I coordinate" i “ time- translation Killing vector and .12b; 

' ' \ t is proper time at radial infinity / 

The gravitational field is characterized by three fundamental variables 
which are functions of and t: 


r [cm] 


M„[gl 


("radius") 3 (l/2«) x (circumference around center of star); (3a) 

( "total mass inside radius r" — Including contributions from rest\ 
mass, nuclear binding energy, internal energy, and gravity 

(3c) 


2 

1 ’ o' 1 - ("8>^®vltatlonal potential") = c^ In [■« 

•(' 

^)l 

Lsec^J 

r 



Energy transport through the star is characterized by three quantities, 
each of which is determined by measurements using standard, physical rods and 
clocks in the mean local rest frame of the baryons at radius r: 


L 

r 





("local luminosity' 


..) - ( 


non-neutrino energy being transported across' 
the sphere at radius r, per unit time 

(4.) 




^ X /-.ame as above, but for neutrino energy\ 

neu r no urn nos y \ ^ . produced in thermonuclear reaction j ; 
rom nuc ear urn ng j cycles which change the abundances X, / 


(" 


neutrino luminosity not 
from nuclear burning 




y cycles which changf 

(4b) 

I same as above, but for neutrino v 

energy produced by processes \ . (4c) 

which, in time- average, do not I 
change the abundances / 


The complete stellar structure and evolution are characterized by the 
functions p(M^,t), T(M^,t), X^(M^,t), r(M^,t), Mj.^(M^,t), 0(M^,t), L^(Mj.,t), 
b"''(M^,t), and L°^M^,t). 
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111. AUXILIARY VARIABLES 

a) Thermodynamic, Nuclear Burning, and Opacity Variables 


The following auxiliary variables are algebraic functions of the funda- 
mental variables; and like the fundamental variables they are determined by 
measurements using standard, physical rods and clocks In the mean local rest 
frame of the baryons 

P(p,T,X ) g (total pressure); (Sa) 

L cm -J 



/binding energy of nuclei, per unit 
I rest mass, relative to hydrogen 



2 


(5b) 


where m^ Is the mass of atomic species 1 in its ground state, m^ Is the mass 
of atomic hydrogen, A^ is the number of baryons in atomic species i, and c is 
the speed of light; 


n(p,T,x^) ^ 


("specific internal energy") 


.( 


total mass-energy of a sample of ] 
stellar material, in energy units/ ^ \ 

(total rest mass of the sample) ' i' 


- c 


(5c) 


^ T V ^ P el /density of total non-gravltational\ „/ 2 „/ 2y /^.v 

Pt(p.T,X^) ^ ^ mass-energy, in mass units j ^ +Il/c ); (5d) 


(p,T,X^) ^ (opacity) = (Rosscland mean opacity); 


(5e) 


h 



I 


'rate, per unit rest mass, at which nuclear^ 
, burning creates non-neutrino energy i 


/ - \ erg per unit rest mass, at which nuclear^ 

VPf » g sec^ ® ^ burning creates neutrino energy ^ 


[i^] 


^late, per unit rest mass, at which non-nuclear-\ 

burning processes [processes with no change ) ; (5h) 

' in create neutrino energy ' 


«i(p>T,X^) j^sec“^j 


( rate at which the abundance X. of species i\ 
changes due to nuclear burning ) * 


b) Relativistic Correction Functions 


The above auxiliary variables (except B and p^) are all familiar from 
the Newtonian theory of stellar interiors. In the relativistic theory it is 
useful to introduce the following additional ancillary variables, each of 
which is dimensionless and is unity in the Newtonian limit 

R H ("redshift correction factor") = exp(o/c^); 

Vs ("volume correction factor") s (l - 2GM^^c^r)""^^^ ; 


^ /"gravitational-acceleration' 


correction factor' 


ation\ 

" j 


+4itr^p/c^ 
tr ' 


e 5 ("energy correction factor") s 1 + (II-B)/c = P^/p} 
Us ("enthalpy correction factor") s 1 + (n - B + p/p)/c . 


In terms of these variables, the general relativistic metric for spacetime 
inside and around the star is 


ds^ = -R^c^dt^ 4 ir^dr^ + r^(dQ^ 4 sin'^Od^^). 


2 > 
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c) Mtxtn g -Length Variables 

The Newtonian mixing- length theory of convective energy transport is 
readily generalized to general relativity. One need only introduce the 
local proper reference frame of an observer at test at radius r, and in that 
reference frame analyze, in a manner identical to Newtonian theory, the 
buoyant forces on convective cells and the heat exchange between convective 
cells and their surroundings. The auxiliary variables that enter into such 
an analysis, patterned after Paczynskl's (1969) Newtonian variant, are 


[ cm * 1 /local acceleration^ 
sec^J ^ V gravity 


\ ^”r 

' r 


Hp [cm] 


= (pressure scale height) = (P/pg)i^^^ ; 


1^ [cm] s (mixing length [normally chosen equal to Hp]); 


u) s (optical thickness of one scale height) = k p 


s (specific heat at constant 


7o s (coefficient of heat exchange) = — 



(8a) 


(8b) 

to Hp]); 

(8c) 

) = K p fj.; 

(8d) 

(bTl\ P lb 

J (Oe) 

'P,Xi 

1 + u)^/3 

(8f) 


b> 


where a s ac/k is the Stephan- Boltzmann constant; and 
O i£L£\ 


(8g) 


In terras of these auxiliary variables, the basic algebraic equations of the 


G 


mixing- length theory are these: (i) An equation which defines the "radiative 

gradient" 


^ /value that (ofnl/ciMy)^ (d/nP/c^M^)^ dfnT/dfnP 

rad ^ would have if the material Were non-convective 




the equation for follows from equations (3.11-3) and (3.11-7a) of Thorne 

(1966, 1967) by straightforward change of notation: 

K L_ P 


3 l./n e\ 

^ V ■ >? ) 


(9a) 


(li) The usual equation for the "adiabatic gradient" 

/ainT\ 




rad 


(9b) 


'entropy,X^ 


where Pg is the adiabatic index of the second kind, (ili) A set of four 
coupled algebraic equations which determine the energy flux carried by con- 
vection F , the mean velocity of a convective cell ("turbulent velocity") 
conv^ 

v^, the gradient associated with a convective cell v', and the actual gradient 
averaged over all convective cells and over the medium through which they 
move, y: 


’'conv-J.pJ, (’rad-'')' 

(10a) 


(10b) 


(10c) 


(v-y')/(y' -V^j) = 7 ^v^ . (lOd) 

Equaclons (lOb)-(lOd) have Identically the same form as in Newtonian theory 
because their derivation in the proper reference frame of a static observer 
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is Identical to that of Newtonian theory. Equation (lOa) also has standard 
Newtonian fonn. I* follows from equation (9a) with rewritten as 
Uitr and from the analogous equation for the actual gradient 

V In terms of the radiative flux 


3 


K(4nr^F 


rad 


)P 


GM oT 
r 






Because equations (lO) all have the same form as In Newtonian theory, one 

can use the standard technique [eqs. (22)-(27) of Paczynskl (1969)] to solve 

them for the four unknowns F , v^, 7', and 7. 

conv' t' * 


IV. DIFFERENTIAL EQUATIONS OF STELLAiv STRUCTURE 

There are 6-t-N (where N Is the number of nuclear species) differential 
equations of stellar structure for the 8 + N fundamental variables p, T, 
r, M , C. L , and as functions of M and t. In these differential 

equations acts at fixed t, and S/5t acts at fixed Each equation Is 

a translation of the Indicated combination of equations from Thorne ( 1966b, 

1967). 

The equation for as a proper volume Integral of p; translation of 
equation (3.11-1): 

dr/^M^ = (Unr^pTr)"^. (lla) 

The equation for total mass-energy inside radius r; translation of 
equations (3.11-2) and (3.11-1): 

e/y. (11b) 
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The source equation for the gravitational potential translation of 


equations (3.11-4) and (3.II-1): 


CM 

r 4nr p 


(11c) 


The equation of energy generation; translation of equations (3.11-5), 
(3.11-6), and (3.11-1); 


5S 


€ - t - — + 

nuc ov 


^ dn . P 1 c>c 
R 5t 


r 1 op 
2 R ^ 


(lid) 


The equation for neutrino losses due to nuclear burning; translation 
of equations (3.11-6) and (3.11-1), specialized to nuclear-burning neutrinos 


1 MlJV) 


R 


"5m" 


at e 


nv 


(lie) 


The equation for non-nuclear-burning neutrino losses; translation of 
equations (3.11-6) and 3.11-1), specialized to non-nuclear-burning neutrinos 


1 Ml°V) 


R 


"5?r 


'ov 


(Ilf) 


The equation for changes of nuclear abundances due to nuclear burning; 
translation of equation ("*,8) 


.-1 


dX^/c)t a . 


dig) 


The equation of energy transport; follows directly from the definition 
of 7; translation of the mixing- length-generalized version of equations 
(3.11-7) 


^ In t/SM^ = V D In P/dMp 


If 

^rad ^ ^ad' 

If 

^rad ^ ^ad' 


(llh) 
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The Oppenheimer-Volkof 1. equation of hydrostatic equilibrium; translation 
of equations (3.11-3) and (3.11-1) 

r 4«r 

This equation must be combined with the equation of state P(p,T,X^) and with 
equations (llg,h) for dT/dM^ and to yield 


V. BOUIT'^RY CONDITIONS 

Corresponding to each different derivative with respect to in the 
equations of stellar structure there is a radial boundary condition. The 
obvious boundary conditions at the star's center are 

r a M, » L « l"'' » L°^' « 0 at M - 0 (l2a) 

tr r r r r ' ' 

(translation of [3.38a]). 

We shall denote the surface values of rest mass, total mass^ radius, and 
the total luminosities by 

nv nv ov ov t v 

L«L,L «L ,L «L at surface. (13) 
r' t tr' ' r' r r 

At the surface the star's spacetime geometry (?) must match onto th* external 
Schwarzschild geometry 

ds^ -(l-2GMj./c^r)c^dt^ + ( 1 - 2GMj_/c^r)"^dr^ + r^(dO^ + sin^ 0 d(p^) . (l4) 

Smoothness of the match ("continuity of intrinsic geometry of surface") 
requires that satisfy the surface boundary condition 

c^(l - PCNj/c^R) at = M (surface of star). (l2b) 
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Note that the luminosities as measured far from the star — which we denote 

nv ov nv 

JLt ^ f fit'll not the same as the surface lumlrov .ties L, L , and 

L°'^. Rather, they are the surface luminosities corrected for gravitational* 

redshlf t 

tt/L - - (1 - 2GM^/c^R) . (15) 

In addition to the boundary conditions (l2a,b) one must also impose 
surface boundary conditions on pressure P and temperature T. If moderate 
errors near the surface are allowable, one can Impose the "zero boundary 
conditions" 

P=T-0 at (12c) 

(translation of eqs. [3.38c,d]). If higher accuracy is desired oi e can 
impose the boundary conditions of the relativistic version of the Eddington 
approximation 

L = 4nR^or , kP “ f Mj. = M (12C) 

(translation of eqs. [3.38c ' ,d ']) . For still higher accuracy one can join 
onto a model stellar atmosphere. If the atmosphere is thin compared to the 
stellar radius R, then it can be constructed in the standard Newtonian manner 
using a surface gravity of 

gg = (GM^/r^) r at = M, (16) 

a surface luminosity equal to L, and radial and time coordinates t and t 
related to r and t by 

t r. (r-R)r, t = tR; fl = = ( 1 - 2GM^/c^R) (l7) 


11 


» 


( 

^ If the atmosphere is not thin compared to R, one can construct it using the 

i formalism of general relativistic radiative transfer theory, which is reviewed 

in §2.6 of Novikov and Thorne (1973). In the case of a thin atmosphere all. 
spectral features as observed by a distant observer are redshifted relative 
to their rest wavelengths by 

= (1-2GM^/c^R)" 2 _i, (13) 


VI. SOME USEFUL RELATIONS 


In this section we list several useful relations among the stellar- 
interior variables. 

The sura of the fractional abundances must be unity at all times; 
and consequently, the sum of their rates of change must vanish: 


E X. = 1, E a. = 0. (19) 

i ^ i ^ 

The total rate of energy release by nuclear burning muft equal the rate 
of change of nuclear binding energy 


c + e 
nuc nv 


1 bB 
Si 'Si 


m^^c 

i "’h^i 


a. 


( 20 ) 


see equations (5b) and (llg). 

The rate of change of the total mass-energy inside radius r, as measured 
by an observer there, must be equal to the rate at w’.-.-'.uh matter carries mass- 
energy inward minus the rate at which luminosity carries it outward: 
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This mass-energy conservation law requires some discussion: (i) The time 

derivatives here are taken at fixed radius r, whereas all previous time 
derivatives were taken at fixed rest mass the two types of time deriva- 

« 

tives are related by 

(s/dt)^ = (a/at) + (aM^/at)^ (a/an^^ . ( 82 ) 

(ii) The operator R is derivative with respect to the proper time 

of an observer who sits at rest at radius rj see equation (?). (iii) 

5? is the locally measured rate at which rest mass flows inward 

across radius r; and (^M is the rate of inflow of rest mass 
plus enthalpy in mass units. Enthalpy appears in the* conservation law 
rather than energy (v rather than e) for the same reason as it appears in 
the Bernoulli equation: in moving matter, pressure (the difference betv;een 

pV and pc) transports energy 

(energy flux) = (pressure) x (velocity). 

(iv) ( l/c^) (L^ + + L°^) is the locally measured rate at which mass-energy is 

transported outward by neutrinos, photons, and diffusive heat flow, (v) Since 
is the total mass-energy inside radius r, one would have expected the 
left side of equation (21) to read R ^ (^M^^/St)^ — i.e., one would nave 
expected the y to be absent. The presence of V suggests to me that 
relativity theor ^ sts such as Misner, Thorne, and Wheeler (1975) should n ot 
have given the name "total mass-energy inside radius r" to . Rather, 
the quantity 

Slitr = (l-r"^) = (c^/G)r |^1 - ( 1 - 2GMj.^/c^r) j 

(23) 

r) p 

— * h '/c'r in Newtonian limit 

tr ^ tr ' 

r.TJ TTY of the 
, I i -X>i' 


15 


should have been Identified as total mass-energy Inside radius r because 


It satisfies 


(^i " ■ 7 ' 21 ') 

Out of deference to established convention I suggest that people retain the 
name "total mass-energy Inside radius r" for but keep In mind that It 

Is a misnomer. 

The equation of mass-energy conservation (21) can be derived from the 
equations of stellar structure by first deriving the relation 




(24) 


where ^/br acts at fixed time t, and by then Invoking the boundary conditions 


n V ov 

L = L = L = M = =0 

r r r r tr 


at r = 0. 


A derivation of equation (24) proceeds as follows; (l) By combining equations 
(lla,d,e,f), ( 20 ), and (Gd) derive the relation 


^ [7 ('•r - 1.7 . = -U„2p„ [(If) - ilf- (If) 

L r r-J 


( 11 ) Use equations (22) to convert from time derivatives at fixed to time 
derivatives at fixed r; and then use equations (lla,c,l) and (6a,d,e) to obtain 


^ [7 (>-r' [(If)^. - (sfjj 




(25a) 


l)^ 



I 


(lit) Use equation (lla) to derive the relation 




(Va) + 



(25b) 


(iv) Use equations (6a,b,c,d,e) and (lla,b,c) to derive the relation 


ib/br)(VR) = 4rt(G/c^)rpVr^B 


and then use equations (lla^b) and (6b) to ibtaln 





(25c) 


(v) Finally, combine equations (25a, b,c) and (6d,e) to obtain equation (214-). 

Note that the equation of mass-energy conservation (21), when evaluated far 
outside the star, just says that If the rest mass of the star is held fixed 
then its total mass-energy decreases at a rate given by the photon and neutrino 
mass-energy losses 

dH^/dt = - {l/c^)U. + (26) 


VII. SUMMARY 

Coordinates M^, t for the stellar interior are defined in equations (?a,b). 
The star's structure is described by 8+N (where N is the number of nuclear 
species) fundamental variables p, T, r, 1*^ » and , which 

are functions of M^ and t, and which are defined in equations (l), (3), and 
(4). These 8 + N variables satisfy the 8 + N differential equations of stellar 
structure (ll), subject to the radial boundary conditions (12). The differen- 
tial equations (ll) contain a number of auxiliary variables, which are 
algebraic functions of the fundamental variables, and which are defined in 


15 


equations (5)-(l0). Quantities which characterize the surface of the star, 
its external gravitational field, and the radiation which leaves the star 
are described by equations (13)-(16). Several useful relations among the 
stellar variables are given in equations (l9)-(2l) and (26). 

This version of the equations of stellar structure and evolution reduces 
to the standard Newtonian version when one sets the following relativistic 

p 

correction factors to unity: R, 7r, <#, e, V in the interior; (l - 2GM^/c R) 

and (l-2GM^/c^r) at the surface and in the exterior. 
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